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Knowledge of a body 's gravitational field and of its anom- 
alies is a boundary constraint for any physical model of its 
structure and the behavior of its small deformations. Such know- 
ledge has now become indispensable to other disciplines contrib- 
uting to the study of the body: seismology, magnetism, 

topography, etc. 

The techniques that have been developed up to now and the 
methods used to determine the parameters representing spatial 
variations, or even space-time variations in the field, have 
produced models of varying accuracy and resolution, not only 
for the earth, but for other telluric planets and the moon. 

They permit comparative structural studies, whose impact on the 
understanding of the geophysical phenomena characteristic of 
our planet has been great. However, even though these models 
appear sophisticated to the general public, they are still 
insufficient to correctly interpret phenomena that affect our 
planet and certain others in the solar system. Probably one 
of the most important of such phenomena is the manner in which 
the variations in gravity, among other geophysical quantities, 
are linked to the plate tectonics and earthquakes. In order 

♦Numbers in the margin indicate pagination in the foreign text. 
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to understand and build models of such phenomena, an understanding 
that is at once general, detailed, and precise of the anomalies 
in the various existing fields, particularly in the gravitational 
field, is necessary. 

The only way to achieve a general knowledge and, later, 
to monitor the time variations in the active zones, involves 
use of outer space. For terrestrial applications, the presence 
of the atmosphere necessitates launching satellites whose 
perturbations are then aralyzed or which carry instruments that 
directly measure varii<"-^ons in the field at an altitude at which 
the effects to be detected are already very weak. Today, the 
space techniques are therefore complemented by the use of grav- 
ity measurements on the ground (or even in tae sea) and by 
altimetric measurements of the form of the privileged equipotential 
surface forming the geoid in its oceanic sections. In the 
oceanic regions, the needs of geophysical studies could, if 
necessary, be met by the joint usage of the overall information 
learned through analyzing satellite orbits and of altimetric 
data. In contrast, the knowledge of anomalies at medium scale, 
and even more so at high resolution, is imprecise in continental 
regions (with the exception of western Europe and part of North 
America) because of poor coverage (often non-existent) and of 
low accuracy inthe measurements of gravity. No improvement can 
be hoped for in this situation, at least as far as the models 
available to the international community are concerned. 

We therefore must look for another method besides alti- 
metry that has equivalent resolution and accuracy, i.e. 
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approximately 1 mgal (= 10 m/sec ) of accuracy for mean 
values in 100 x 100 )an regions. 

There essentially exist two space techniques to achieve 
these goals, satellite- to-satellite tracking and gradiometry. 
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Satellite-to-satelllte tracking has already undergone 
study In France (Project Dlabolo) , Europe (Project Slalom) , and 
the United States (Project Gravsat) . This method, based on 
the analysis of the relative motion of at least two satellites, 
Is certainly very complicated and costly. Its feasibility has 
not yet been proved. On the other hand, gradlometry Is con- 
ceptually a simpler method (only one satellite) , but It has 
come up against difficulties that seem to put the necessary 
accuracy outside our grasp. However, the success of space accel- 
erometry end ONERA's recent research have both Indicated a way 
to overcome these obstacles. The first studies and experimen- 
tation carried out by ONERA and ONES are very encouraging and 
call for a much deeper Investigation of the mission and system. 

In the United States, the development of space gradlometry Is 
linked to the use of ground technologies which carry with them 
heavy penalties when used In outer space. On the contrary, we 
are proposing here to start off directly from a preexisting 
space technology: the CACTUS accelerometer built by ONERA and 

launched by CNES on board the CASTOR satellite In 1975. The 
situation Is therefore favorable for meeting current, attainable 
scientific objectives. 

1. Scientific Objectives 

1.1 Knowledge of the Earth's Field 


A more and more detailed global knowledge of the earth's 
gravitational field Is necessary for several purposes: 

— Improvement In reconstructing and extrapolating the 
trajectories of artificial satellites for practical purposes 
(navigation, precise positioning by quasl-geometrlc methods — 
dynamic or short arcs) and scientific ones, such as analysis of 
residual perturbing effects on the trajectories for the study, 
e.g., earth and oceanic tides; determination of polar motion 
by methods using radio measurements (Doppler or laser) ; 
reduction of altlmetrlc measurements. 
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— Use of altlmetric data , not only in geodesy , where this 
information is employed to define a mean reference source close 
to the geoid (and in certain cases equal to it ) , but also in 
space oceanographies , provided that this mean surface can be 
exactly distinguished from the geoid and its temporal variations 
studied. 

In geophysics , at various scales. Figure 1 (GRIM 3 Geoid) 
shows a recently-made chart of the earth's overall field as 
determined by analyzing the perturbations of geodetic satellites 
(observed between 1960 and 1970) and by comparing them with 
surface information: gravitational anomalies and altlmetric 

heights of the geoid. This model, the third of its type in 
Europe, is the product of Franco-German cooperation. It is 
representative of the long and medium wavelengths in the grav- 
itational field (X/2 ^ 500 )cm) , as the geoid represented here 
shows, and corresponds to the harmonics of potential up to an 
order and degree of thirty- six. The modulations in the geoid 
(or associated gravitational anomalies) in relation to a 
reference ellipse express the divergence between the figure and 
the earth if the earth was in hydrostatic equilibrium. They 
are directly linlced to lateral heterogeneities in internal 
density. The goal of interpreting the anomalies is to determine 

the nature of the internal constraints associated with the 

variations in density, for example, to see whether the earth is /J_ 

sufficiently rigid at the depths studied to behave in an elastic 
fashion or whether, if it has a more plastic behavior, a 
dynamic support is necessary (by means of convective motion) . 

The large-scale undulations in the geoid (X/2 > 1500 ]on) 
do not show any evident correlation with the topography, 
whether continental or oceanic. This suggests that they 
reflect anomalies in deep density, located a few hundred )cilo- 
meters down in the mantle, or even more. Such anomalies would 
be related to convection cells within the mantle, proposed as 
the motive forces of the observed movement in the tectonic plates. 
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Despite numerous experiments and various mathematical models, 
the phen(xnenon is still poorly understood and suffers frcxn 
errors which continue to affect the gravitational field in many 
regions. This is illustrated in figure 2 , which depicts the 
differences between the geoid model cited above and another 
model (American) that is practically equivalent. The divergence 
is especially great in the continental zones and is mainly due 
to the errors tainting the gravitational anomalies used (or to 
the absence of surface information.'). 

For shorter spatial scales, the relationship between 
gravity and topography, especially the marine topography of the 
great oceanic structures, is remarkable. A striking correlation 
exists between the geoid 's medium wavelength undulations and 
the tectonics of the ocean floor. Each type of structure, sub- 
duction zones, dorsal regions, fracture zones, and continental 
shelves, has its own particular signature that has been well 
detected by altimetric techniques. This perfectly expresses 
the fact that isotasy works at these wavelengths. It is 
obviously desirable to study such compensation at this scale, 
and with at least the same accuracy, on the continents, so as 
to correlate medium undulations with measurements of the crust's 
horizontal and vertical deplacements. This is where space 
geodesy has a great role to play in the observation of active 
zones. It can be predicted — but is impossible to fix a 
date! — that the precise determination of ^e gravitational 
anomalies in these regions, the measurement of the variations, 
combined with the continuous measurement of horizontal and 
vertical displacements, will lead to an understanding of the 
mechanisms that trigger earthquakes and, perhaps, contribute 
to predicting them. 

The indirect techniques used up to now have practically 
reached their limits, even when eexobined with surface information. 
This is because the problems in the coverage and accuracy of 
gravimetric measurements cannot be overcome in places where 
altimetry is indispensable for reasons of accessibility (for 


example, the Himalayas) , politics, and so forth. In addition, 
even if the earth were covered by oceans, their mean surface, 
as determined by altimetry, would not coincide with the geoid, 
and these are the very differences that interest oceanography. 

We therefore have to come up with an alternative technique. 

There seem to be only two possibilities at present, satellite* 
to-satellite tracting (SST) and gradiometry. 

The SST technique consists of measuring the radial velocity 
between two satellites in one of the following configurations: 
one satellite high (geosynchronous in scmie cases) and the other 
low (HL) ; two satellites in low orbit (LL configuration) and 
close to one another (5'J to 300 km) . HL was investigated using 
GEOS 3 and ATS 6, but confirmation of the technique's validity 
only occurred once (GEOS 3 was too high, at approximatley 800 
km) . LL was attempted during the Apollo-Soyuz rendezvous, but 
the test was not conclusive because of the two vessels' ccmiplex 
spurious motion and ionospheric disturbances. In each case, 
the measurement is interpreted as a difference in gravitational 
potential between the two objects. The better technique seems 
to be the LL configuration, and several projects in this 
direction have been designed; 1) Project Slalom (formerly 
Diabolo) investigated first by the GRGS alone and then in 
collaboration with the SFB78 and MBB within the framework of 
the ESA. This project is based on laser techniques. It uses 
two Spacelab subsatellites, but the many difficulties linked to 
the space shuttle seem to be insoluble . . . and the mission was 
very limited, both in space and time. 2) Project Gravsat-A 
in the US^ consisting of two enormous satellites (6 m long 
cylinders 1.50 m in diameter, three-axis stabilized, and con- 
taining several tons of propellant so as to ensure a one year's 
life in orbit between 160 and 190 km) and measuring radial 
velocity b> radioelectric technique (approximately 110 GHz) . 

The feasibility of the American project (studied by NASA and 
the DOD in collaboration with (illegible word)) is far fran 
proven. Even though its estimated cost, which is already very 
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high, does not seem to be an insuperable obstacle (for the DC 
the experiment cannot take place before 1987 at the earliest. 
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fig ire 3 

Minimum Instrument Accuracy ior SST/LL and Gradicxnetry Techniques 
as a Function of Orbital Altitude and for Resolutions of 
1* X !• and 1.5" x 1.5" 

Key: a) Average Orbital Altitude (Km.) 

b) SST/LL 

c) Gradicmsetry 



The gradi<^etrlc technique is to measure one or several 
ccmiponents of the tensor of second derivatives of gravitational 
potential in the axes linked to the satellite (there are five 
independent components) . It has never been the object of an 
experiment in outer space. Under study since the end '^f the 
sixties in the USA, gradicxnetry has led to the design and some- 
times the construction and testing (oit the ground) of prototypes. 
There presently exist four American gradicxnetry projects which, 

in the presence of gravity, are reported to be capaL .* of achieving 
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a resolution of 0.1 EU (1 Ebtv&s Unit > 10 sec , corresponding to 
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a field of 1 mgal « 10 cm-sec in 10 )cm) . However, all of these 
projects suffer frcxn unresolved problems either with the instruments 
or the environment, and none is yet ready for a space flight. At 
best, the most recent of them (H.J. Paik's cryogenic gradiometer) 
might be tested in the space shuttle after 1985 and carried by 
satellites towards the end of the decade. 

In view of the scientific objectives we havw enumerated, a 
resolution of X/2 in the 110 to 165 km range (1® x 1® and 1.5® x 1.5® 
respectively) could be fixed. An accuracy ii: reconstructing 
gravitational anomalies at this scale on the order of one milligal 
could likewise be decided upon. We can then determine the altitude 
of the satellite or satellites and the accuracy of the instruments 
necessary to attain these objectives. Figure 3 summarizes the con- 
clusions of an analytic study founded on an extrapolation to short 
wavelengths of the orders of magnitude of ancmialies in potential for 
both techniques, SST/LL and gradiometry. 

As far as gradi(xnetry is concerned, the reader can therefore 
bear in mind the ccxnprcxnise given in the following table: 
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Table 1 

Precision of Gradiometric Measurement (EU) as a Function of 
Altitude for the Two Resolutions Given 


ti tude 

RSsoluti ori\^ 

200 

250 

300 

350 

400 

r X 1“ 

0.10 

0.025 

0.006 

0.002 

0.0004 

l?5xl?5 ' 

0.40 

0.20 

0.075 

0.030 

0.015 


We shall see in the following chapters how we can meet /12 

such conditions. 

1.2 In Planetology 

The gravitational field of earth-lika planets remains the 
information of choice for constructing any model, even an 
elementary one, of internal structure in the absence of direct 
(seismic) measurements. Comparing both the planets' fields with 
that of the earth and the structures of large and medixim scale 
anomalies with topographical features and surface tectonics 
is a way to investigate these bodies' origins and thermal and 
physio-chemical history. 

Mars, outside of the regions whose latitude is less than 
55® S, is certainly the best known from this point of view, thanks 
to the Mariner 9 and Viking 1 and 2 missions and to the orbits 
with very low periastrons from which we have benefited. The 
resolution of the detected anomalies is on the order of 500 km. 
Considering the size of the planet, this is similar to global 
models of the earth. Gradiometry is probably the only technique 
that could appreciably improve resolution, if not coverage, 
given the constraints in possible future missions. 


Venus offers more Interesting possibilities. This steins 
from the fact that the only sophisticated mission engaged in 
this type of study, the Venus Pioneer, did not in fact resolve 
the anomalies sufficiently except in the neighborhood of the 
periastron of the orbiter (whose latitude was almost constant) , 
i.e. for a zone at 40® latitude. The only coming mission which 
might carry enough apparatus to answer questions concerning 
tectonics and isostasy on Venus is VOIR. However, VOIR, as it 
is envisioned today, requires that experiments be of long 
duration due to the observational constraints arising from the 
planet's slow rotation. Whether a gradiometer can function on 
a vehicle of the type planned and alongside the other instruments 
remains to be proved. 

Of course, we can contemplate missions to Mercury, but 
this planet is outside Europe's capabilities in view of launch 
constraints and the small ntunber of scientific instruments it 
would be possible to carry on board. At present, the United 
States is not planning any mission to this planet. 

The earth-like body which is scientifically the most inter- 
esting from the point of view of origin, evolution, and, why not, 
future development is our natural satellite. The ESA is studying 
a very complete mission around the moon called Polo. It has 
taken up the main themes of the American project LPO, which was 
abandoned at the end of the seventies for budgetary reasons. 

One of its objectives is to produce a complete cartography of 
surface physio-chemistry, magnetism, topographical relief, and 
gravitational anomalies. For this last subject, the planning 
is to set up a relay satellite in order to use the SST/HL tech- 
nique to determine the gravitational field of the hidden side. 
This approach increases the mission's cost and complexity 
considerably. It is clear that however one wishes to continue 
studying this mission, i.e. within the ESA or in direct bilateral 
cooperation with NASA (which is in the process of redefining 
a new, ambitious lunar mission for the end of the decade), the 


elegant and imposing solution is to place a gradiometer on 
board the orbiter. The fact that the experiment takes place 
between 50 and 100 km above the lunar surface means that even 
here a precision of 1 EU, at most 0.1 £U, is satisfactory, 
given the expected anomalies. 

Finally, and this can be considered an extreme case that 
is not without interest, a gradiometer could "simply" serve to 
measure the mass of a body such as an asteroid. For example, 

3 

an asteroid with a diameter of 200 km and a density of 3 g/cm 
would produce a gravitational gradient of 8 EU on a gradiometer 
flying by at an altitude of 500 km, a gradient of 26 EU at 300 
km, and so forth. In the first case again, the signal would 
be 27 EU if the asteroid had a diameter of 300 km. Such measure- 
ment would be fundamental during an asteroid flyby mission 
given the large uncertainty (a factor of two) concerning the 
density of the materials that make up asteroids. 

1.3 Other Applications of the Technique 

The basis of a gradiometer ' s operation lies in determining 
the relative motion of two (or more) small test masses placed 
inside cavities, without any material contact with the fixed 
cavities on board the satellite. Because of this fact, the 
test masses do not undergo any of the nongravitational surface 
forces .^.ting the satellite (atmospheric braking and 
radiation pressure, for example) . The restoring forces, of 
electrostatic origin, exerted on the test masses to keep them near 
the center of the cavities therefore make it possible to 
determine these surface forces. In the mission envisaged 
(circular polar satellite at around 250 )an or lower) , the total 
density of the atmosphere can thus be determined very precisely. 
Even though knowledge of this value has made great progress 
since the beginning of the space age, its long-term evolution 
over different solar cycles is still not well enough known. 

Exact knowledge of it would be very useful for the kind of 
satellite reentry envisaged in several programs. 



g. Instrumentation; ONERa's Concept of the Gradiometer 


The purpose of the gradiometer is to measure the tensor of the 
second derivative of the potential of gravitational attraction. The 
technique is presently developed, particularly in the United States, 
for high resolution determinations of gravitational potential. 

In the forseeable applications with satellites flying by at 

300 km from the earth* s surface and 100 km from the moon, for example, 

the difference between the gravitational attracti'^n at two points 
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1 m meter away from each other is around 10 m/sec . This difference 

.... -9 -2 

is equivalent to 1000 EU. (1 EU = 1 Eotvos Unit = 10 sec .) 

A gradient of this magnitude suggests that a gradiometer based on a 
system of ultrasensitive accelerometers like ONERA's CACTUS [Ultra- 
sensitive Three-Axis Capacitive Accelerometric Sensor] accelerometer 
could be built. 

So as to pinpoint the present state of the technology and make 
the report easier, we shall review the description and use of this 
type of accelerometer before taking up the gradiometer ' s operating 
principles . 

2 . 1 The CACTUS and SUPERCACTUS Accelerometer [1 to 7] 


2.1.1 History 

ONERA specially developed these accelerometers for space 
applications. Thus, CACTUS made up the payload of the French 
satellite CASTOR that the CNES developed and launched on May 17, 1975. 
The data acquired during its 45 month orbital life were analyzed 
by the Grasse Center for Geodynamic and Astronomic Studies and 
Research (CERGA) . The first scientific studies concerned atmospheric 
temperature and density based on the braking experienced by the 
satellite between 270 and 400 km altitude. 
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While cactus's expected accuracy was 10 m/sec , statistical /16 
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processing of the data allowed an accuracy of 10 m/sec to be 
achieved. 

Thus, it was in particular possible to; 

— follow the annual variations in solar radiation pressure; 

— measure and construct a model of the earth's albedo; 

— measure the acceleration due to the earth's infrared 
radiation. 

Finally, the analyses documented well the thrust the satellite 
felt from its own thermal radiation, whose anisotropy was due to 
variations in surface temperature. When the satellite was exposed 
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to the sun, this thrust showed up as an acceleration of 3 x 10 
m/sec . It went down to zero within a few minutes after the satel- 
lite entered the earth's shadow. 

SUPERCACTUS is a new, higher resolution version, a preliminary 
study of which was made on the request of the European Space Agency. 

Use of this accelerometer should make a satellite for measuring all 
the various components of terrestrial radiation (Project BIRAMIS) 
possible. 

2.1.2 Description and Model of the Accelerometer 

The operating principle of CACTUS accelerometers consists of 
measuring the force developed by triaxial electrostatic suspension 
to maintain a spherical inertial mass -- a ball — in the center of 
a cage attached to a satellite. 

Figure 4 is a schematic representation of this servocontrol ; 

f designates the force coming from the satellite that is 
exerted on the ball of mass m. (? is a binding force.) 
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fg is 
applied to 


the resultant of the nongravitatlonal exterior forces 
the satellite of mass M - m. 



Figure 4 

Schematic Diagram of How a Sensitive Axis Functions 

Key: a) Detection Electrode c) Telemetry 

b) Action Electrode d) Servocontrol Electronics 

Figure 5 specifies the notation used: In the XYZ frame of /17 

reference the vectors and define the positions of the satel- 
lite's center of mass and of the center B of the accelerometric 
ball. A servomechanism keeps the ball at the accelerometer's center 
0. Consequently, in steady state operation the positional divergence 
in oS is negligible. The ball's position is then defined in relation 
to the satellite's frame of reference by the constant 

vector Using and Tt to represent the gravitational accel- 

eration at G and B, the equations of motion in the fixed reference 
system for these two points is written: 
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Figure 5 


Notation: 

X Y Z: fixed reference system 

X y z ; satellite's frame of reference 
s s s 

G: satellite's center of mass 
0: center of accelerometer 
B: center of ball 



or for 
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To operate the accelerometer and determine I’g, it is very 
desirable to "center” the satellite so that G and 0 occupy the same 
position. Since in steady state operation the ball's center B is 
equal to 0, the three points G, 0, and B all coincide. Equation (1) 
is then reduced to: 



m, K 


such that knowledge of f makes it possible to directly determine the 
resultant of the nongravitational forces acting on the satellite. 

In the general case and when several accelerometers are used, /18 
the points G and B are distinct. To take account of acceleration 
due to the satellite's angular motion, it is then convenient to 
express equation (1) in the form it takes when projected onto the 
satellite's frame of reference x^y^z„. if the time derivative in 
this frame of reference is designated by a and n used to 

represent the instantaneous rotational vector, it occurs that: 
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( 2 ) 
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More generally, if n accelerometers are arranged about the 
satellite, one obtains for each of them: 


^E^i, + ^ 4- JTI v>i, 

a 

M _ £ «i 


•llNCSKSBi') 

* 5^i ■ 


Since in steady state operation each is kept in position at the 
center 0^ of the corresponding accelerometer, it occurs that: 


^3;, r and ^ _ 


such that each accelerometer's output signal is finally 

expressed as: 
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2.2 Principle of the Gradlometer 


2.2.1 Gradiometry Using Two Accelerators 

If the value to be measured is the difference ( 
between the gravitiational accelration at G and 0^, equation (3) 
shows that the accelerometer's output signal is not sufficient. 

It must be freed of the satellite's acceleration 


fe - 


as well as of the acceleration due to its own angular motion. For 
that it is natural to turn first to the difference in the measure- 
ments made by two accelerometers. Taking equation (3) into account/ 
one can directly obtain; 

£ - fi. s r -V - 

The differential measurement eliminates the effect of the 

external force I*-, but it still remains to evaluate the inertial 

E 

acceleration 


o 


For 0/1 Oa, U « 0.5 m, an angular velocity of 2 X lO"^ 

rad/sec or an angular acceleration of 4 x 10 rad/sec along 

the directions perpendicular to 0,6_ would give rise to an acceler- 
• 6 2 * * 

ation of 2 X 10 m/sec . Such a perturbation is equivalent to 
the action of a gravity gradient of 4000 EU. 


of: 


In other words, to obtain 0.1 EU accuracy requires estimation 


— SXx , with an accuracy of , or, for the example 

* -8 
chosen, to the nearest 2.5 x 10 rad/sec; and 

O 10 2 

- j:i, to che nearest 10~ rad/sec^. 

In practice, this level of accuracy is not attainable with the /20 
apparatus used for reconstructing satellite attitude. The gradio- 
meter alone must therefore nullify the effects due to the satellite's 
angular motion . This requires at least four accelerometers. 

2.2.2 Gradiometry Using Four Accelerometers 

The basic principle consists of defining a particular arrange- 
ment such that inertial acceleration can be eliminated when 
processing the accelerometric measurements. Figure 6 illustrates 
the arrangement envisioned. The centers 0^, O^f 0^, of the 
four accelerometers are arranged at the vertices of a unit square. 

The center 0 of the square does not necessarily coincide with the 
satellite's center of mass G. 

In what follows, the vector equations are written directly in 
the instrument's frame of reference xyz. JCI and St describe 
this reference system's angular motion relative to the fixed one. It 
is necessary at this stage to make the gravity gradients explicit. By 
definition the corresponding tensor (T) is symmetric and its trace 
is zero (the derived field of a potential) . 


(4) 



’’'■4*1 



A 



T44 


with 


ix<x 




s. O 



Figure 6 

Gradiometry Using Four Accelerometers: Arrangement Adopted 

X, y, z: the gradiometer platform's frame of reference 

0: center of the platform 

0^: center of accelerometer i 

G: satellite's center of mass 

platform's angular motion 

In all, therefore, (T) only includes five independent elements. 

In considering the acceleration at the center 0 of the square 
constituted by the gradiometer, it happens that; 



and it is possible to distinguish two terms in equation (3) of the 
measurement produced by each accelerometer: 



JL • r ♦ rj, 


f is an acceleration, acting in common mode on all four 
accelerometers : 


(7) Pa ^ ^ 4 . — ^ ^ 4- II- ^ SO 4 . 


7^ is jm acceleration directly proportional to the sepa- 
ration . 




The square arrangement results in: 


OO3 m - OOi, 


which makes: 


t; 


= - r. 


(9) 


and 
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Then, considering the body of measurements supplied by the 
four accelerometers: 


I 


( 10 ) 



i varying from 1 to 4, one obtains: 


( 11 ) 



/22 


Equation (10) ma)ces it possible to operate a 
together as an accelerometer. Equations (11) are 
usable for the gradiometer and lead to a total of 
which are easily written out by ta)cing account of 


syst^ working 
the only ones 
six scalar equations, 
( 8 ): 



24 



one obtain;3: 

/ 


ORIGINAL . 

OF POOR QUALITY 


( 12 ) 


"Htx. - '* 

Tx^ • ^ ■*' 

Tx^ - ■*• SlxSl^- j 

T*x **■ ^->t.<Sl^ - .2^ - 

* t ^ 4 ^ ^ - ^Tt ’♦’ 



This system includes 11 unknowns (fl, n, (T)) and is only linear 

if the direction of the instantaneous vector of rotation is known . 

It is therefore desirable to impose this direction in a way that 

->■ 

minimizes the errors due to the uncertainties in nandi'j. For that, 
it has proved to be particularly interesting to orient the satel- 
lite's angular velocity n along the z-axis perpendicular to the 
plane of the square Thus; 


•9 ^ 




T 


which supposes in addition that: 


(14) 


Ja » [ o 
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System 12 now has only seven unknowns; 

o 

- and on the one hand, and 

— the five independent elements of tensor (T) on the other. 


After eliminating £ 1 ^ and it should therefore be possible 

to obtain four independent equations between the gravity gradients 
from the six equations in (12) . 
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However, conditions (13) and (14) are never perfectly satisfied, 

o A 

and, in particular, the residual values of Ax. are directly 

expressed as errors in the gradients. For accelerations and 

• “lO 2 “10 2 

of 10 rad/sec , the error in Txz (or Tyz) is 10 sec , 

or 0.1 EU. This leads to restricting the measurements to the com- 
ponents located in the gradiometer ' s xy-plane. There are only 


three tensor elements active in this plane (Txx, Txy, and Tyy) . 

After eliminating and the four equations to retain in system 
(12) finally maJce it possible to obtain two relationships containing 


these three elements. Using 
residual values of and 
calculations : 


An* and to represent the 

, one obtains, after all the 



The effects due to the angular acceleration SI are completely /24 

eliminated. The residual velocities AA« and only intervene 

in the form of squares or products such that an accuracy of 0.1 EU 
-10 -2 

= 10 sec corresponds to the specifications: 


tsi. 


LSI. 
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The attitude control system has to be able to accommodate 
this constraint. 

2.2.3 Implementing the System and Anticipated Accuracy 

When the velocity deviations and are zero and the 

square perfectly regular, the accuracy is only limited by: 

— the satellite's contribution to the gravity gradients Txx, 

Tyy, and Txy, and 

— the accuracy of the accelerometers themsexves. 

The satellite's contribution to the gradients may attain 10 EU, 
but presents the peculiarity of being practically constant in the 
gradiometer ' s xyz reference system (which is connected to the 
satellite). The satellite's gravitational attraction has to be 
evaluated in flight. While not attributable to the operating scheme 
adopted for the gradiometer, it should be taken into account when 
designing the satellite and analyzing its mission. 

As for the accelerometers' accuracy, the error in measurement 
obtained along a sensitive axis mainly results from three terms. For 
the x-axis of accelerometer i, the measurement is F ./m. and the 

XI 1 

error is expressed as: 


(16) 







I 




is a bias resulting from spurious forces inside the accel- /25 

erometer and whose anticipated value in the SUPERCACTUS is lO”^^ 

2 

m/sec . For 1 = 0.5 m, this value corresponds to an uncertainty on 
the order of 0.2 EU. Thus, not only is this effect industinguishable 
from the satellite's contribution to the gradients, but, in addition, 
its value is much smaller. 
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represents the thermal noise emitted by the electronic 

servocontrol circuitry. This noise is mostly found in the highest 

frequencies of the accelerometer's passband and, as a result, can be 

easily filtered out. Thus, for SUPERCACTUS, the use of a low-pass 

-12 2 

filter with a 15 sec time constant reduces bxi to 10 ra/sec , or 
for 1 = 0.5 m, a noise equivalent to approximately 2 x 10~^ EU. 

This value represents the instrument's thermodynamic limit when the 
integration time is 15 sec. It is advisable to remember that in 
spite of the design adopted, the inertial acceleration arising from 
the satellite's angular motion in reality causes perturbations that 
are a priori greater than this theoretical limit. 


AS Txi 

X — expresses the error corresponding to the relative 

uncertainty in the accelerometer's scale factor S. Since the 

-4 2 

accelerometer's range is too small (less than 10 m/sec ) to allow 
it to function on the ground, pre-liftoff determination of its 
sensitivity S can only be accomplished indirectly. Under these 
conditions, it is not possible to anticipate an accuracy of better 

_3 

than 10 for AS/S. Thus, in the presence of acceleration F . /m. on 

—6 2 ^ 
the order of 10 m/sec , the value obtained should be; 




-9 2 

“ 10 m/sec , or an uncertainty in the 
gravitational gradients of about 2 EU (for 1 = 0.5 m). 


Even though it is desirable to minimize the value of the 
acceleration F ./m., equation (3) shows that the limit to consider 
is fixed by | * 1 value is precisely on the order of 

10~® m/sec^. It should be noted in addition that, among the accel- 
erations endured by the instrviment, corresponds to the amount of 

atmospheric braking exerted on the satellite, which varies greatly 
with altitude. 


Since the values of the acceleration are only susceptible to /26 

being increased, one must try to reduce the error produced by the 
uncertainty AS/S in the scale factor. Now, gradiometry relies only 
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on differential measurements. Thus, when the same systematic error 

in sensitivity taints all the accelerometric measurements, this error 

transfers directly onto the gradient measurements (for the same 

reason as an uncertainty in the length 1 of the sides of the squares 

does) . Attainment of accuracies better than 1 EU (for accelerations 
— 6 2 

of 10 m/sec ) therefore requires that the accelerometers ' scale 
factors be equalized in flight . Such an operation could be performed 
by a servocontrol system like the one represented in figure 7. 

A calibration device generates a sinusoidal acceleration r 
at pulsatance a> in the accelerometers' passband. This acceleration, 
which can be obtained by rotating a small, unbalanced wheel connected 
to the gradiometer at constant velocity, for example, is applied 
in common mode to both accelerometers A^^ and A 2 whose sensitivities 
and S 2 are to be equalized. Synchronously demodulating A^^'s and 
A 2 *s output signal at a pulsatance w makes it possible to extract the 

difference (r - r_) between the measurements obtained for r. The 

JL * 

sensitivity S 2 is then adjusted so as to equal S^ by introducing a 
modification AS 2 proportional to the integral of 
last value remains constant when S^ = S 2 * 

It is difficult to define beforehand the limit of the accuracy 

-4 -5 

in this feedback control. Values on the order of 10 to 10 seem 
reasonable. 

In conclusion, the evaluation of the gradiometer ' s overall 

accuracy should be the occasion for a serious parametric study of 

the instrument, the satellite, and the mission. However, the 

envisaged operating principles and the orders of magnitude presented 

indicate that a resolution of better than 0.1 EU in the presence of 

_6 2 

perturbing accelerations at least equal to 10 m/sec is possible. 
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Figure 7 

Servocontrol for Equalizing Scale Factors S, and S 
of Accelerometers and A2 

a) Generation of a sinusoidal acceleration 

b) External Acceleration 

c) Synchronous demodulation 

d) Modification of S~: 







3. The Terrestrial Mission 
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3.1 Technical Analysis of the Mission [8] 

3.1.1 Principle of the Mission 

To obtain a good representation of the earth's potential with 
the aid of sufficiently exact measurements of a certain number of 
components of the tensor of the earth's gravitational gradient. 

The measurement of these components is performed by making 
three-axis differential acceleration measurements with SUPERCACTUS 
type accelerometers arranged according to a certain geometric con- 
figuration. 

The basic configuration used is one with accelerometers dis- 
tributed about a square (figure 8) . 

Obviously, the earth mission is much trickier than the other 
planetological ones accessible to the envisaged gradiometer because 
of the earth's atmosphere and the precision sought (beyond the 
36th harmonic) . 

3.1.2 Main Constraints 


— Global coverage of the earth in small enough steps. 

— No saturation of any microaccelerometer's output (which would 
entail the total loss of the information) . 

— The spurious acceleration terms must-, be identified with 

a precision compatible with the resolution sought after, in 
conformity with table 1. 
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Figure 8 

Geometry and Arrangement of an Instrumental Group 
Relative to a Local Orbital Frame of Reference 

Key: a) velocity 

b) normal to the orbit 

c) center of the earth 

3.1.3 Discussion of Principal Factors 

Orbit : Almost total coverage of the earth requires an orbit of 

high inclination and therefore one that is practically polar. In 
addition, the opportunities for sending up secondary passenger 
transports (Viking philosophy) into heliosynchronous orbits during 
launchings of multimission platforms of the SPOT variety support 
the choice of such an orbit as the parking orbit. Taking account of 
the cost in mass of major orbital corrections, a final orbit that 
is accessible at the end of a reasonable period of time and with a 
reasonable increase in speed is called for. 
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We will therefore assume in what follows that we are working /31 

in a heliosynchronous or quasi-heliosynchronous orbit. This 
results in coverage of more than 99% of the earth's surface (99.36% 
for 96.5“) with two parameters to determine: altitude and local 

time of node. 

Local time of node : The preliminary analysis of the problems 

of measurement and both attitude and orbital control led to recom- 
mending a measurement system composed of floating gyros and suggested 
the use of ionic engines to partially compensate for atmospheric 
drag. The attitude and orbital control system will therefore consume 
significant amounts of energy. The obvious necessity to have a very 
low A/m (this point is developed in detail below) is the basis for /32 

rejecting solar generators that turn or are too large. 

To produce enough power with a generator stuck to the vehicle's 
skin, a major main frame member must point towards the sun. In order 
to minimize A/m in the direction of the aerodynamic wind, this 
maximum main frame member must be perpendicular to the normal of the 
orbit. The sun's directon must be as close as possible to the orbit. 

In conclusion, the local time of the ascending node has to be 6:00 or 
18:00. The choice depends on the primary orbit dictated by the 
principal passenger (minimization of drift time) . 

Altitude . This parameter affects: 

— The continuous components of acceleration (aerodynamic drag, 
continuous vehicle rotation, first te’im of the potential) 
which threaten to saturate the accelerometer. 

— The sensitivity required for isoresolution gradient measure- 
ment in anomalies (table 1 and figure 3) . 

— The quality of coverage in terms of step, time, and redun- 
dancy in a short mission (six months or one year) . Table 2 
presents the principal parameters related to the quality of 
coverage. 


33 


The useful lifespan in terms of maintenance within the range 
of altitudes in which measurement is possible. (The useful 
lifespan is sharply lower than the lifespan in the traditional 
sense of the term.) 


Table 2 

Coverage as a Function of Chosen Orbit 


Mean Altitude 

(ka) 

Inclination* 

(o) 

Orbital Period* 
(mn) 

Coverage 
Resolution (kn) 

** 

in n days 

200 

96.408 

88.395 

220 


24 




84 

- 

48 




32 

- 

95 




20 

- 

190 

220 

96.477 

88.798 

280 


16 . 




90 

- 

46 




16 

- 

230 

. 240 

96.546 

89.202 

330 


18 




80 

- 

67 




10 

- 

472 

260 

96.616 

89.607 

300 


19 




MO 

- 

37 




35 


112 

280 

96.687 

90.012 

64 


78 




6 


390 

300 

96.758 

90.418 

400 


12 




130 

- 

35 




22 


210 


3 2 

*Calculations made with R = 6370 km, GM = 39,8601 km /sec , 
C 2 Q = -1.082628 x 10-3, eScifitricity held at 0.001. 

**According to the case, cert5in regions might in fact be better 
covered than the table indicates. The figures therefore give the 
niaximum distance between two traces. 
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The first analyses carried out showed that above 300 km measure- 
ment accuracy is limited by the accuracy in reconstructing the angular 
velocities that induce the inertial accelerations * Qr + n^r, where 


us 

-«y 


MX 

My 

o 

-MX 

MX 

0 


that the gradiometer cannot distinguish from the acceleration gradients 
of the earth's potential. 

In addition, taking account of the predictable gradiometer dynamic /34 
range, the vertical channel threatens to saturate should the vehicle 
rotate at the orbital velocity. 

A preliminary analysis of this point is found in [9]. 

The simplified system of equations giving the measurements' 
accuracy is: 


A ( T - T ) < Aty + 4« .Am 

XX yy m 

AT < 2 Ay > 2m. Am 

xy m 


where : 


— AVm is the gradiometer ' s error in measurement expressed as 

a function of the measurement field y i Ay - cy . 

o m o 

— M is the residual angular velocity (maximum model). 

— Aw is the gyro unit's error in reconstructing the angular 
velocity. 
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As a first approximation, the values adopted for u and Lu are 
-4 >4 

respectively 5 x 10 and 2 x 10 deg/sec. Starting from there, 
the choice of measurement field and dynamic range 1/c establishes 
the measurement accuracy in Eotvbs units. This leads to the maximum 
altitude as determined with the aid of figure 3. 

For the problem to have a solution, the accelerations acting on 
each accelerometer and each channel (aerodynamic drag, solar radiation 
pressure, gravity gradient of the spherical earth, and inertial forces 
connected with rotation relative to the Inertial reference system) 
must be less than y^. 

If the vehicle rotates about the normal to the orbit, the dom- /35 
inant accelerations are: 

— Along the local vertical YY; 

Vertical gravitational gradient 

2 

Rotation about the normal t aiz 

— Along the tangent XX: 

Horizontal gravitational gradient 

2 

Rotation about the normal t uz 

Atmospheric drag 

— Along the normal ZZ: 

Horizontal gravitational gradient 

Solar radiation pressure (negligible) 

The first evaluations reveal feasibility problems when using a 
theoretically realistic value for the gradiometer's dynamic range, 
i.e. e = lO”^. 
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This suggests the use of an ionic engine to at least partially 
compensate for atmospheric drag> this XX term being the most 
critical. The use of such a propulsive system would in any case be 
necessary to c(»npensate for erosion of the orbit at this altitude 
(even with the lowA/mhoped for), since the range of allowed altitudes 
for correct measurements is small. 

Furthermore, the level of thrust is limited to 10 mM 
because of the ionic engine's high electrical consumption (375 W 
for 10 mN) . 

If T is the drag and F the engine's thrust, it is advisable 

to ensure that 1 f - t!/M < y in order to have measurements and that 
o 

F - t > 0 along an orbit. 

Reconstructing attitude and orbit ; The accuracy in reconstruc- /36 
ting the attitude is expected to be on the order of a half-degree 
(effect less than one-tenth EU) , and in this situation the control 
accuracy is not important. 

As for the orbit, the principal constraint concerns reconstruc- 
ting the altitude with an accuracy on the order of 10 to 100 uteters 
so as not to introduce too great an error as altitude decreases. 

3.1.4 Summary of Preliminary Conclusions 

The measurements themselves ; The accuracy and dynamic range 
of the measurements have to be studied more in depth. This particular- 
ly concerns; 

— the final accuracy in control and reconstruction of angular 
velocities as a functionof the passband; 

— the exact value of the gradicxneter ' s dynamic range 1/c 
(equalization of scale factors, especially) as a function 
of the passband; 
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— th« spectral nature of the desired information so as to judge 
its spe'iral separation trom different perturbations. 

T ho orbit ! Preferable orbit: 6:00 or 18:00 heliosynchrcnous 

oroit with an altitude of between 200 and 300 km. Reconstruction of 
vertical position: to the nearest 10 to 100 m. 

The vehicle should have: 

— a minimum A/m perpendicular to the aerodynamic wind, 

— a maximum surface area perpendicular to the normal to the 
orbit; 

— an almost constant attitude relative to the local orbital 

frame of reference: orbital control rotation at orbital 

pulsatance; 

— very low angular velocities (flight control with high 
damping coefficient, elevated inertia, reduction of perturb- 
ations) and the lowest possible pulsatance, and measure them 
with precision; 




— a capability to make major orbital changes; /37 

— an ability to change attitude for doing calibrations or 
making complementary measurements; 


— provisions for providing the instrument itself with a high 
degree of mechanical and thermal stability (between the four 
accelerometers and between the instrument and the vehicle's 
frame of reference) . 


38 



ORIGINAL page: ;S 





Figure 9 

The vehicle's external shape is lenticular (or discoid). It is 
centered on a cylindrical collar of the standard 1200 diameter joining 
the SPOT multimission platform to the Ariane (exactly like on the 
Viking) and has a height of 0.5 to 0.55 m. 

The type of control system that governs the vehicle's flight 

relative to the local orbital refe’*ence svstem allows for a minimum 

2 

main frame member (less than 1 m ) "acing the aerodynamic wind and 

2 

a maximum one (greater than 3 m ) facing the sun. 

Considering the Ariane 3's performance in an orbit of the SPOT 
type (822 km), i.e. 3400 kg, and the predicted maximum mass of 2000 kg 
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for a satellite using the SPOT multimission platform, we can hope for 

a mass of around 1000 kg and therefore a ballistic coefficient 

A/m of 10 (Remember that theDSB— = 5.6 x 10~^.) 

m 

The mass allowance over and above the preliminary tally of mass 

requirements is enough to plan for an annular ballast of heavy 

material around the edge to achieve a high level of inertia, greater 

2 

than several hundred m -kg. This margin of unallocated mass also /38 

has to cover the lower braces connected to the central tube and the 
1200 diameter adapter. A more detailed analysis has to be made of 
how the central tube can be arranged to allow for the cohabitation 
of the instrument and the propulsive system. 

3.2.2 Mechanical and Thermal Architecture 


The essential points of the structural design concern: 

— the braces connecting the peripheral ring and the 
central tube; 

— the connection between the instrument and the central tube 
(dimensional and directional stability) . 

On the subject of heat distribution, the vehicle will have the 
side covered with solar cells continually facing the sun while the 
opposite side will be in the shade. It is therefore advisable to 
provide a high degree of thermal coupling between the two sides. 

In addition, the radial and circumferential gradients have to be 
minimized (or at least their variation has to be) because they 
directly affect the bias due to the satellite. 

3.2.3 Attitude and Orbital Control System 

Normal Mode [9] 


The recommended attitude measurement system includes a gyro unit 
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derived (by simplification) from that of the SPOT and two photocell 
strip star trackers. The attitude control systCTi should be based 
on the use of flaps controlled by incremental motors. A set of four 
flaps that are symmetric in relation to the xz and xy planes and 
located in a plane parallel to the y-axis will provide three-axis 
attitude control, particularly since the other perturbing torques 
are negligible in comparison with the aerodynamic torque. 

Orbital control (compensation for drag in order to measure and 
maintain altitude) will be accomplished with the help of an RIT type 
ionic engine (manufacturer; MBB) whose thrust can be adjusted 
between 5 and 10 mN before takeoff. 


Its specifications are; 
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Thrust; 

I ; 
sp 

Fuel Supply Mass; 

Engine mass; 

Electrical Consumption; 


between 5 and 10 mN 
31,000 NS/kg 
11 kg 
4.3 kg 

230 W for 5 mN, 375 W for 10 mN 


Future studies will develop a strategy for using this propulsive 
system that takes its high electrical consumption into account. 


Attitude acquisition ; The satellite will be placed in an 800 km 
circular orbit with a random initial attitude and an initial velocity 
that might reach 12 deg/sec about each axis. Achieving the nominal 
attitude will therefore require an attitude measurement 
system. 


Attitude acquisition will take place in the following phases; 

— reduction of velocity components; 

— solar pointing with the aid of a fine solar sensor; 

— Commencement of slow spin about the pointed axis; acquisition 
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of stars by the star tracker and telemetric processing, with 
evaluation of the satellite's attitude using ephemerides; 

''I 

— transition to operational mode by using the star trackers. 

During these different phases, the attitude control Systran will 
have to produce torques along the three axes with minimal fuel 
consumption. Considering the relatively low inertias, the torques 
will not have to be too great. 

A cold gas jet system therefore seems appropriate. This is /40 

especially true since it would allow attitudinal maneuvers in the 
800 km orbit for possible instrument calibrations and could induce 
the spacecraft to spin or alter its uctitude in preparation for orbit 
changing impulses. 

Mass and power tally sheet (outside of the system for changing 
orbits) ; 



Acquisition 

Normal Mode | 


mass 

power 

mass 

power 

Attitude Measurement 

21 kg 


45 W 

21 kg 

45 W 


Orbital and Attitude 

( - 



1 X 16.3 

375 W 


Control 

( 6 kg 

{ - 


6 N 

4 kg 

0 W 


Propellants 

3 kg 


— 

22 kg 

— 


Total 

73 kg 

51 W 

73 kg 

420 W 
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Orbit acquisition ; The vehicle will have to move from a circular /41 
orbit to one with a different plane and altitude. With the exception 
of proposals involving aerodynamic braking, -\t least two impulsions 
will be necessary. Furthermore, the satellite will have to pass 
through an intermediate drift orbit because a direct correction would 
be too costly in terms of mass. 

For an 800 km orbit, the drift relative to heliosynchronism is 
given in 110] ; 

^ (hours per day) = -7 x 10~^Ai (degrees) + 1.4 x 10~^Aip (km) 

Cthe orbit being an elliptical one with an apogee of 800 km) . 

The first correction simultaneously alters the semimajor axis 
and the inclination. With the perigee lowered to 300 km, Arp equals 
-500 km, and the associated negative drift term has a value of 

^ = -7 X 10 ^ hours/day, 

“(rp) 

One possible strategy is to transfer into a 300 x 800 orbit by 
changing inclination, drift for 270 days, and then make the orbit a 
circular 300 km one while maintaining the inclination at the heliosyn- 
chronous value. 

— First correction; Arp = 500 km Ai = +1.38® 

AV = 225 m/sec 

— Drift for 270 days in a 300 x 800 km orbit, i = 100®. 

— Second correction; Ara = -500 km Ai = -3.4® /42 

AV = 460 m/sec 

The two corrections taken together correspond to a propulsive 
mass equal to 21% of the initial mass. 
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It will be necessary to make a deeper analysis of this strategy 
and to choose between a system of two solid-fuel engines and a liquid 
fuel system. 

3.2.4 Generation of Power 


The angle A between the sun and the normal to the orbit varies 
between +17® and -30®, according to the season. Thus, eclipses occur 
in one part of the year (A > 17®) . This part of the year is centered 
on one of the solstices. Which one depends on the local time chosen 
for the node. 

The electrical power available from the generating system and 
originating in the cells glued to one whole side is about 270 to 314 W, 
according to the season, with a diameter of 2 meters. Depending on 
the operating cycle adopted for the ionic engine, it is possible that 
it will be necessary to add a crown of cells around the edge. A 
crown 25 cm deep will make the available power Increase to between 
420 and 490 W. 

3.2.5 Other Subsystems 

The other subsystems do not pose any particular problem provided 
that a more precise evaluation of the aerials necessary for the 
TM-TC-LOC system that takes account of the desired localization 
accuracy (10 to 100 m) is made. 

3 . 3 Ground Segment /43 

Thw problems linked to the ground segment are of two orders: 

— recovery of the telemetric data stored on board; 

— decametric localization of the satellite to permit reduction 
of the measurements. 
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The amount of data to transmit per orbit (about 6 Mbits) does 
not present any obvious problems for the ground stations currently 
available or for the use of bubble memories. 

For reconstructing altitude, the measurement of drag, the use of 
a model with a large enough potential (see TOPEX study) and of a 
few laser stations should solve the problem. 

3.4 Development Plan 

A very preliminary evaluation ended in the following results: 
Feasibility Studies (Phase A) 

Duration of about 1 1/2 years, end of 1981 to mid 1983. 

Essentially concerns: 

— Scientific mission: precisely define the scientific missions, 
simulate the instrument's operation in the adopted config- 
urations, evaluate the system's overall accuracy. 

— Instrument: gradiometer feasibility study, limit of resolution, 
specifications of inflight calibration mechanisms. 

— Mission and vehicle: detailed analysis of mission, choice 
of orbit, orbit acquisition, identification of possible 
configurations, general architecture of satellite, atti- 
tudinal movements and deformations, corrections, stabili- 
zation. In principle, the subsystems should not pose any /44 

feasibility problems since we are building on equipment 
that is either already in existance or presently under 
development (SPOT gyros, Symphonie cold gas system, photo- 
cell strip star trackers of the ERS 1 type) . 

— Elaboration of the development plan and evaluation of its 
cost. 
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Detailed Design (Phase B) 


Duration of about two years, mid 1983 to mid 1985. 

During this phase, the plan is to conduct complementary technical 
studies of SUPERCACTUS and to create functional mockups of the gradio- 
meter and the inflight calibration mechanisms. As far as the system 
is concerned: detailed study of its functions (attitude control, 

orbital control, thermal control, electrical generation, data 
acquisition and transmission, structure, mechanical architecture) , 
detailed design of the satellite, ground stations and data processing, 
system and subsystem specification, organization of the construction 
phase, precise determination of costs. 

Development and Construction (Phases C and D) 

Duration of about 2 1/2 years, mid 1985 to end of 1987. 

Customary model fabrication and qualification phase, particularly 
qualification of SUPERCACTUS and the complete platform, equipped 
with calibration mechanisms. 

The studies in phases A and B have to determine which preliminary 
models to build ahead of the complete qualification and flight models. 

Key Dates /45 

The principal stages might be the following (See chart below) : 

— End of 1981: Decision to launch phase A. 

-- Mid 1982: Synthesis of the first feasibility study results 

(particularly the first simulations) and decision to pursue 
phase A. 

— Mid 1983: Review of phase A and decision to launch phase B. 
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3.5.1 Measurement Reduction 

Among the problems connected to measurement reduction, one merits 
particular attention because of its possible impact on the satellite's 
design. The test masses are subjected to both internal gravitation 
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forces, linked to the satellite (distribution of masses inside the 
satellite) , and external ones (earth, moon, sun) . We shall always 
be analyzing the resultant of these forces. However, the purpose of 
the mission concerns analyzing the external forces. To separate the 
effects, one must work in different conditions; reduction of the 
data obtained at very different altitudes or reduction of data 
collected with permutations in the role of the satellite's axes. One 
will therefore consider: 

— the advantage to a temporary intermediate orbit, for example 
at 800 km; 

— the advantage to and the necessity of being able to rotate 
the satellite about itself in order to exchange the axes (as 
was done with CASTOR) . 

These problems must be studied with care. 

3.5.2 Calculation of Field Anomalies 


Any determination of parameters linked to the geopotential at 
the earth's surface, such as geoid altitude or gravitational anomalies, 
based on observations made on board a satellite suffers from what 
the geodesists call the "problem of downward prolongation" . This is 
simply the problem of the attenuation of the signal caused by an 
anomaly with increasing altitude. It is illustrated by table 2, where 
we have tabulated the variance of gravitational anomalies (extra- 
polated beyond an order and degree of 36) for blocks of wavelengths 
(i.e. between two given degrees and as a function of altitude. 

We could say, for example, that the signal is considerably /47 

amplified in the 73 to 180 block when one descends from 250 km to 
0 km, and so is the error in the detected signal, unfortrnately . 

This is the first difficulty (almost nonexistent up to £ = 36, as all 
the present overall models prove a posteriori ) which does not yet 
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have a satisfactory theoretical solution. However, the smaller 
the error in the measured components is, the less troublescxne it will 
be. This consideration might justify a (redundant) system to measure 
the six tensor components. In fact, only fine covariance analyses or 
very realistic simulations can show the limits in reconstructing the 
parameters we are interested in. In any case, a first observation 
made by all the theoreticians is that by measuring the components 
of V T Cwhere T is the perturbing potential) , the propagation of 
errors is less for the calculation of gravitational anomalies 
{'v- I |vt| 1) than for T itself, or in other words for tracing the 
geoid. This difficulty has already been pointed out. 

Table 2 

o (Ag; -*■ ^ 2 ^ mgals, for Different Altitudes 


^■■v.,,Spectral Extent 

; f 2 

Altitude 

(km) 

2 18 

19 - 36 

37 72 

73 180 

400 

8.85 

2.35 

0.69 

0.06 

250 

10.82 

4.25 

1.96 

0.40 

0 

15.97 

12.44 

14.09 

17.54 


The reconstruction process, for example of gravitational anomalies, 
raises problems of another nature. The use of spherical harmonics 
is impracticable because this is not a regionalized representation of 
the field. Point masses or surface densities could be used, but it 
is even more preferable to use the generalized Sto)ces equation, which 
presents an even more "pointed" regionalization and malces it possible 
to directly obtain the gravitational anomalies near the subsatellite 
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trace. In theory, the relationship connecting the (infinite) vector 
of values of g to the (finite) vector of measurenient q Is discretized /48 
In the form: 


where the operator S Implicitly contains Plzettl's extension of the 
core of the Sto)ces operator. Mathonatlcally , ll)ce all Inverse prob 
lems, this problem Is called "Improperly stated", and a possible 
solution consists in carrying out a TDchonov regularization by 
looking for solutions g such that 


where D is the covariance matrix of measurement errors, Cgg is the 
(infinite) covariance matrix of the field g (estimated in steps or 
extrapolated from what one already knows of the selfcorrelation 
function of g at the present wavelengths) , and a is a regularization 
factor. The solution has to be judged by its stability in relation 
to "reasonable" variations in a or perturbations in the matrix Cgg. 

In any case, it must be recognized that the practical resolution 
of such a problem requires sophisticated computing apparatus and 
might involve dividing up the steps of the calculation within Europe, 
even at the final level. This is similar to what we have already been 
doing for the last decade with the C. Reigber group in Munich. 

4. Conclusion /49 


Significant progress in the knowledge of the gravitational 
field of the earth and other telluric bodies necessarily passes by 
way of the development of new measurement techniques, either satellite- 
to-satellite tracking or gradiometry. 
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Gradiometry, in essence the simpler of the two, now seems within 
our ability to perform, and the first analyses, which are performed 
and briefly presented in this document, show the possiblity of 
occupying an original niche in this field. 

The preliminary design is based on a gradiometer with four ONERA 
microaccelercaneters. It calls for an approximately 1,000 kg lens- 
shaped satellite stabilized in the local orbital frame of reference. 
Its orbit would be a 6:00-18:00 heliosynchronous one ranging between 
200 and 300 km in altitude. 

This mission to a large extent uses the French space experience 
acquired in the D SB (Castor) and SPOT progr 2 uns. 

In contrast, it appears that we have insufficient knowledge of 
the spectral nature of gravitation measurements and their possible 
separation from perturbations, the intrinsic accuracy of the 
gradiometer ' s measurements (equalization of scale factors), the 
accuracy in controlling and reconstructing the vehicle's angular 
velocity, and the methods of dealing with the data (model creation) . 
These various points do not, at the present stage, bring the first 
conclusions into question, but they do necessitate cc»nplementary 
studies to o.nfirm the project's feasibility and the ultimate accuracy 
attainable in the measurement of ancxnalies in the gravitational field. 
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